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The  magnetohydrodynamic  (MHD)  equations  that  describe  the  motion  of  a 
tenuous  plasma  in  a  uniform  magnetic  field  are  mathematically  analogous  to 
the  hydrodynamic  equations  describing  two-dimensional  flow  of  an  incompress¬ 
ible  fluid.  As  a  consequence,  one  can  construct  a  laboratory  analogue  of 
MHD  phenomena  involved  in  the  structuring  of  high-aititude  plasma  clouds  on 
scales  that  lead  to  radiowave  scintillation.  A  simple  laboiatory  tool  for 
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visualizing  the  relevant  dynamics  is  the  Hele-Shaw  cell,  which  consists  of  two 
pieces  of  plate  glass  between  which  appropriate  fluids  are  constrained  to 
undergo  two-dimensional  flow.  Two  fluids  of  different  density  but  equal 
kinematic  viscosity  are  placed  in  a  Hele-Shaw  cell  with  a  uniform  boundary 
between  them,  and  the  cell  is  oriented  vertically  with  gravity  acting  across 
the  boundary.  Within  a  few  seconds,  the  boundary  becomes  structured  under  the 
influence  of  the  Ray leigh-Tay lor  instability,  which  then  may  be  studied  well 
into  its  nonlinear  regime. 

In  the  dynamic  analogue  that  ensues,  the  fluid  density  and  velocity 
represent  respectively  the  plasma  density  and  convective  velocity  (in  the  E  x 
B  regime),  the  fluid  stream  function  represents  the  electrostatic  potential, 
and  the  gravitationally  driven  free-fall  velocity  in  the  cell  represents  the 
gravitationally  assisted  neutral-wind  velocity  in  the  ionosphere.  In  this 
report,  the  mathematical  underpinnings  of  the  analogue  are  reviewed;  they  are 
found  quite  firm  for  scales  in  excess  of  about  100  meters  in  the  nighttime 
equatorial  F  layer,  subject  to  one  untested  assumption. 

Under  this  contract,  a  Hele-Shaw  cell  capable  ofi^roviding  repeatable 
experimental  conditions  necessary  for  systematic  stiray\  of  Rayleigh-Taylor 
structuring  was  built  and  several  exploratory  experimental  runs  were  made. 
Difficulties  were  encountered  with  a  densitometer  used  for  data  preparation, 
but  useable  results  were  obtained.  The  cell  and  laboratory  procedures  pro¬ 
duced  sufficient  photographic  signal-to-noise  ratios  to  demonstrate  that 
definitive  studies  of  the  spatial  spectrum  of  Rayleigh-Taylor  structure  are 
possible  in  the  spectral  regime  for  which  the  presently  developed  mathe¬ 
matical  analogue  is  valid. 

All  spectra  obtained  show  a  clear  break  between  a  white-noise  regime  and 
a  power-law  tail,  interpreted  as  an  outer  scale.  All  experimental  runs 
analyzed  reveal  an  increase  in  the  outer  scale  near  the  original  fluid  inter¬ 
face  during  their  tens-of-seconds  duration.  Photographic  records  suggest 
rather  nonuniform  spatial  statistics  as  a  function  of  distance  from  the  origi¬ 
nal  interface.  This  report  presents  results  on  evolution  of  the  outer  scale 
near  the  original  interface,  suggests  means  for  probing  its  control,  and  poses 
other  spat ial-stat ist ical  questions  that  could  yield  to  Hele-Shaw  analysis. 
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I.  INTRODUCTION 


For  about  a  decade,  it  has  been  known  that  there  is  a  direct  mathematical  analogy 
between  the  motion  and  structuring  of  a  variable-density  plasma  embedded  in  a  magnetic 
field,  in  the  low  col  1 i s ion-f requency  regime,  and  the  motion  and  structuring  of  a  variable- 
density,  incompressible  fluid  embedded  in  a  resistive  medium  (Thomson,  1970).  Comparison  of 
the  structures  that  develop  across  magnetic  field  lines  in  high-altitude  barium-ion  clouds, 
under  the  influence  of  the  gradient-drift  plasma  instability,  with  those  that  develop  in 
fluid  released  in  a  Hele-Shaw  cell,  under  the  influence  of  the  Ray leigh-Tay lor  instability, 
lends  credence  to  the  analogy. 

Thomson’s  analogue  and  the  experimental  evidence  for  its  validity  offer  a  relatively 
inexpensive  way  to  investigate  some  aspects  of  high-altitude  magnetoplasma  structures  that 
can  disrupt  performance  of  radars  and  communication  systems  by  means  of  radiowave  scintil¬ 
lation.  This  report  summarizes  first  attempts  to  exploit  the  analogue  to  establish  certain 
characteristics  of  the  spatial  spectrum  of  such  structures. 

Work  on  the  contract  has  shown  that  the  Hele-Shaw  technique  is  quite  a  useful  tool  for 
characterizing  the  spectrum  at  scales  larger  than  about  100  m  in  the  topside  ionospheric  F 
layer.  The  one-dimensional  spectrum  found  is  well  characterized  by  a  cutoff  power  law,  and 
the  Hele-Shaw  technique  reveals  dynamically  the  evolution  of  its  outer-scale  cutoff. 

The  work  reported  herein  included  a  review  of  Thomson's  theoretical  development  and 
early  Hele-Shaw  experiments  and,  more  recently,  construction  of  an  improved  cell  and  per¬ 
formance  of  several  experimental  runs  in  it.  Theoretical  considerations  are  presented  in 
Section  II.  Section  III  provides  a  description  of  the  new  cell,  together  with  the  labor¬ 
atory  and  data-recording  techniques  employed  for  experimental  runs  in  it.  Procedures  for 
data  analysis,  including  a  nagging  hardware  problem,  are  described  in  Section  IV.  Results 
and  initial  conclusions  regarding  the  evolution  of  fluid  (and  magnetoplasma)  stuctures  are 
presented  in  Section  V,  along  with  suggestions  for  broader  exploitation  of  the  Hele-Shaw 
technique. 
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II.  THEORETICAL  CONSIDERATIONS 


A.  FIRST  POINT  OF  ANALOGY 

In  order  to  provide  a  sound  basis  for  applying  the  Hele-Shaw  analogue,  a  detailed 
review  of  the  MHD  equations  whose  solution  would  describe  the  evolution  of  plasma  structure 
in  the  ionospere  was  undertaken,  along  with  a  review  of  the  simpler  corresponding  hydro- 
dynamic  equations.  In  an  effort  to  illuminate  all  important  assumptions  and  approximations 
that  might  be  required  in  developing  the  analogue,  the  rather  general  approach  taken  by 
Thomson  (1970)  was  followed.  In  this  section,  we  discuss  the  momentum  equations,  starting 
with  the  following  ion  and  electron  equations  written  in  the  earth  frame: 

dU.  ^  ^  ^  -t 

nm  ■  =  -  VP  +  en[E+(U.  *B)  ]  +  nm.g  -  nm.  v.  (V . -V)  -  m  v.  tj  (1) 

1  dt  l  l  l  l  in  l  e  le  e 


-VP  -  en[E+(U  xB)  ]  +  nm  g  -  nm  v  (P  -V)  +  m  v  ir 
e  e  e  eene  eiee 


where  —  stands  for  the  convective  derivative,  and  where 
d  t 

n  =  plasma  number  density 

m-  =  ion  or  electron  mass 
i  *e 

U-  =  ion  or  electron  velocity 

V  =  neutral -wind  velocity 

P.  =  ion  or  electron  pressure 
i  ,e 

e  =  elementary  charge 

it  =  electric  field 

B  =  magnetic  induction 

g  =  gravitational  acceleration 

v.  =  ion-neutral  collision  frequency 
in 

v  =  electron-neutral  collision  frequency 
en  j 

v.  =  ion-electron  collision  frequency 

w-  y 

and  J  =  ne(U.  -  U  )  -  current  density. 

i  e 

The  beginning  assumptions,  then,  include  local  neutrality  (all  relevant  distances  ^ 

Debye  length),  a  single  ion  species  of  charge  e,  and  infinitely  heavy  neutral  particles 

(which  can  be  accounted  for  in  evaluating  v-  ) 

°  in" 

We  proceed  to  the  first  important  equation  of  the  analogue  by  adding  Eqs.  (1)  and  (2)  to 


VP  +  (0*0)  +  nMg  +  nMy  V  -  n(m.v.  U.+m  v  U  ) 


i  in  i  e  on  e 


where  we  have  defined 

M  =  m .  +  m  =  m 

1  e  l 

tn.U.  +  m  0 

5  =  -1  -1--  e  e  «  0. 

M  l 

P  =  P.  +  P  =  2P  '  2P . 

p  l  e  e  i 

m  v .  +  m  v 

,  i  in  e  en  _ 

and  v  = - — — — —  -  v  . 

n  M  in 


Manipulation  of  the  last  term  in  cn.  ( "1 ) ,  without  approximation,  yields 
dV  -*  mimo  0 

nM  —  =  _  Vr  +  (JxB)  +  nMg  +  nMv  (V-U)  +  Vj  ~ 

dtp  n  M  d 


wher°  v,  =  v  -  v.  :  v 

d  en  in  en 

v  “V 

But  since  m  Svin.  E  while  and  U.  are  comparable  and  v,  is  of  similar  order  to  m 
e  i  e  i  K  d  n 

(although  somewhat  larger),  we  ignore  the  current  term  in  Eq.  (4).  We  then  invoke  Ampere's 

law  in  the  electrostatic  approx (mat  ion,  as  r>''nvs: 

->  ->  1  SE  _► 

V  x  b  =  u  .T  +  — -  ~ ■  ~  u  J 

ccs3to  (a) 


where  -  =  magnetic  oermeabi l i t y  of  free  spac". 

Accordingly,  we  have 

-►  -►  (VxB)xB  _  -VB?  (B « V) B _ A 

•I  ^  B  fy  •  O 

O  0  0  o 


or,  for  a  magnetic  field  up^  fort,  in  its  owr  clitic 


->  ->  V  B 

J  x  B  =  ~ 


where  P  =  - —  is  the  magnetic  pressure.  With  these  further  modest  restrictions,  Eq. 
ib  2  P 

(4)  becomes 
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( 8 ) 


n>‘ 


dU 

dt 


-  VP  +  n?:c  -  n.V v  U 


where  P  =  P  +  P  =  the  total  pressure 
P  m 


and  0  -  g  +  v  v  =  the  acceleration  due  tc  gravity  and  the  neutral  wind.  A 1  terns t i ve 1 v ,  we 
n  •  ’ 

tsv  i pnore  gradients  and  motions  parallel  to  the  magnetic  field  and  write 


n> 


,  dtj.  _ 

dt 


-  V.  P  +  nMC.  -  nMv 

i  l  n  1 


(9) 


ics .  (t;  and  (Q )  are  very  similar  to  the  Nsvier-S tokes  equation  for  an  incompressible 
r  ’  u d  -r  '  C  .  io/1K>  r~_  oc  f„i  lows! 


r 


dU 

dt 


Vp  +  f,g  +  yjV(' 


VU  )  +  It  .  •  u 

0 (in compress ibil 1 ty ) 


(  ^  r  \ 


»”  ?r  ~  ‘"iiit  ve  .  *c  * 

P  =  mass  density 
p  =  hydrostatic  pressure 
11  =  viscosity. 

To  obtain  an  equivalence,  we  must  investigate  the  final  term  in  Eq.  (10)  in  the  Hele-Shaw- 
cell  configuration. 


A  Hele-Shaw  cell  is  simolv  two  pieces  of  plate  glass  separated  by  a  very  small  gap  of 
dimension  d,  between  which  fluid  is  constrained  to  undergo  two-dimensional  flow.  In  the 
analogue,  the  mapne t ic-f ie 1 c  direction  corresponds  to  the  in ter p 1  a te-gap  direction.  Since 
there  is  no  motion  across  the  gap,  the  pressure  and  gravity  terms  in  that  direction  balance, 
and  Ea.  (10)  is  two-Jimensiona 1 ,  as  follows: 


du^ 
dt  ’ 


+ 


nv2u, 


(11) 


where  n  =  —  is  the  kinematic  viscosity.  Further,  the  steady-flow  approximation  is  well 
satisfied  both  in  the  field  and  in  the  laboratory,  so  Eq.  (11)  gives  approximately 


V 


(12) 


For  a  small  gap  dimension,  d, 


we  suppose  that 


.,2 

V  M  ■ 

i  i 


mates 


7  v 

9^  ui 


so  Eq .  ( 12 )  approxi - 


8 


(13) 


±JLL 

•7. 


For  constant  kinematic  viscosity,  which  we  impose  on  our  experimental  fluids,  we  can  write 

(14) 


Oz 


( ■  *4  ■  [> 

( --■£  +  -i’ 

_  '  iP  \ 

\  n  /  n  L1' 

\  IX-  Jy 

.  •  /  *z_ 

and,  from  the  balance  of  forces  across  the  gap,  we  have 
Sp  _ 

57  "  ”  P8„ 


(15) 


_3_ 

3z 


(16) 


For  a  vertically  oriented  cell,  =  0.  If,  also,  we  can  ignore  any  density  variation 
across  the  gap,  then,  the  right-hand  side  of  Eq.  (13)  is  independent  of  z,  and  we  can 
integrate  to  obtain 


^Z-t\ 

V~r2  +  Bz  +  A 


(17) 


where  A  and  B  are  to  be  evaluated  from  the  boundary  conditions  at  z  =  0  and  z  =  d, 
respectively.  Denoting  the  slip  velocity  at  the  boundaries  by  us,  we  obtain  the  following 
cross-gap  velocity  profile: 

?iP 

(18) 


ul  =  us  + 


2n 


(z  -zd) 


which  has  the  value 


■>  d 

uc  ~  us  8 


?iP 


(19) 


at  the  gap  center. 

Under  the  foregoing  conditions,  then, 


9 


(20) 


->  2  ■>  8  .  ->  \ 

'  ul - 2  (uc‘us) 


and,  for  zero  slip  velocity,  Eq.  (11)  becomes 


dt 


(21) 


where  we  have  dropped  the  distinction  between  and  its  maximum  value,  uc.  Rewriting  Eqs. 
(9)  and  (21),  we  have,  for  the  plasma  and  cell  fluid,  respectively, 


dt 


-v 

nM 


+  2  -  v  V  -  0 

i  n  1 


(22a) 


and 


0 


(22b) 


where  n  is 


the  kinematic  viscosity, 


£ 

P* 


Equations  (22a  and  b)  disclose  the  first  point  of  analogy  between  the  Hele-Shaw  fluid 
and  the  magnetoplasma.  They  are  acted  upon  respectively  by  hydrostatic  pressure  and  mag¬ 
netoplasma  pressure  (both  gas-dynamic  and  magnetic);  gravity  drives  the  instability  in  the 
cell  as  can  the  gravitationally  assisted  neutral  wind  in  the  ionosphere,  and  scaled  kine¬ 
matic  viscosity  in  the  cell  plays  the  roll  of  the  ion-neutral  collision  frequency  in  setting 
up  a  drag  force. 


B.  ELECTRODYNAMICS  IN  THE  ANALOGY 

While  Eqs.  (22a  and  b)  do  indeed  suggest  an  analogy  between  Hele-Shaw  fluids  and 
magnetoplasmas  of  interest,  the  very  important  electrodynamics  of  the  former  are  quite 
hidden  in  the  magnetic-pressur*  term  of  Eqs.  (7,  8,  and  22a).  To  bring  them  out,  Thomson 
(1970)  followed  the  procedure  of  considering  electron  and  ion  dynamics  in  detail  to  obtain 
relationships  between  the  density,  J,  of  electric  current  flowing  in  the  plasma  and  the 

-V 

plasma's  convective  velocity,  .  Thus,  if  instead  of  adding  Eqs.  ( 1 )  and  (2 ),  we  subtract  them 
after  dividing  them  respectively  by  nm^  and  nmg ,  we  get 
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Manipulating  various  terms  and  defining 
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but  introducing  no  approximation,  we  obtain 
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On  all  time  scales  of  ^interest,  we  may  ignore  electron  acce  leration,^so  the  left-hand 

dUi  .  dl'i 

side  above  approximates  — —  while  that  of  Eq.  (4)  approximates  ran  — — .  Scaling  and 

dt  dt 

equating  the  right-hand  sides  and  performing  some  manipulations,  while  applying  only  the 

2  2 

innocuous  approximation  mp  <<  nn  ,  we  obtain 
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Equation  (25)  is  valid  even  for  accelerating  ion  flow.  To  proceed  with  our  case  of 
interest,  we  make  the  steady-flow  approximation,  which  is  valid  for  subsonic  flow  and/or 
moderate-to-large  plasma-density  gradients.  With  both  electron  and  ion  acceleration  neg¬ 
lected,  we  may  set  the  left-hand  side  of  Eq.  (24)  to  zero  and  solve  for  the  plasma  velocity, 
U,  as  follows: 


(26) 

where  we  have  defined  the  gyro  frequency 

► 

eB 

“c  =  M"  • 

Inserting  Eq.  (26)  into  Eq.  (25)  and  performing  some  vector  manipulations,  we  obtain 
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where  z  is  a  unit  vector  in  the  direction  of  the  magnetic  field  and  is  the  magnitude  of 
field-aligned  current  density. 

—y 

The  right-hand  side  of  Eq.  (27)  may  be  thought  of  as  an  equivalent  electric  field,  e  , 
driving  the  currents  appearing  on  the  left-hand  side.  Defining  also 
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we  may  write  Eq.  (27)  as 


(rx  +  r2)  J  -  r  J|(  i  +  r3(J  x  z)  «  e  , 


which  may  be  solved  for  the  vector  current  in  terms  of  the  equivalent  driving  field  and  the 
Pederson,  Hall,  and  parallel  conductivities  as  follows: 


at  x  e 

->  c 

J  =  0  e  +  o  -  +  o  e,  z 

p  H  a)  ii  H 

c 


where 
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The  approximation  in  the  conductivity  definitions  holds  for  —  ”  _  <<:  1  and  therefore  is 

2  i  c 

valid  so  long  as  ( v  +  v.  )v.  is  not  >>  u  . 

en  le  in  c 

From  the  continuity  equation  for  charge  density  under  either  the  electro-static  or 
local-neutrality  approximations,  both  of  which  were  invoked  in  developing  the  momentum 
equation,  we  get 

V  •  j  *  0  .  (29) 


Combining  Eqs.  (28)  and  (29)  yields 
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Inserting  the  approximate  forms  for  Fp  and  G  into  the  definition  of  f  and  invoking  —  >?  - — 

>  > 1,  we  find  e  ln 


r  ,  i  VP  /  VP  \  a) 
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Now  employing  Faraday's  law  in  the  magnetostatic  approximation,  we  express  E  in  terms 
of  a  scalar  potential,  <t>  ,  and  write 


:  +  e(^  X  M  +  (e  8  -  Tif )  X  ^ 
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where  we  have  defined 
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With  these  definitions,  the  current  divergence  equation  takes  the  form  of  the  following 
potential  equation: 
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A  series  of  vector  manipulations  without  approximation  leads  to 
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where  we  have  defined 

f  =  ”  vv  «  F  v  V  +  (g-^Ml  • 
enf  e|_n  \s  nM  /  J 

In  the  same  >  >  mg  approximation,  Eq.  (26)  for  the  plasma  velocity  becomes 


(36) 


If  the  current  is  expanded  in  terms  of  the  equivalent  driving  field,  £  ,  and  the  three 

component  conductivities,  °  ,  and  then  Eq.  (34)  is  inserted,  a  rather  complicated  expression 

■+ 

for  U  results.  It  may  be  manipulated,  however,  and  without  introducing  any  new  approxi¬ 
mation,  broken  into  the  following  components  parallel  and  perpendicular  to  the  magnetic 
field: 


In  the  regime  vn  '  “  »  then,  we  have  for  the  perpendicular  velocity 


(37a) 


(37b) 
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In  the  same  regime,  Eq.  (35),  which  expresses  the  fact  that  the  current  is  divergence-free, 
comprises 
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If  we  suppose  the  magnetic  field  strength  to  be  approximately  uniform  in  the  x-y  plane,  we 
may  divide  Eq.  (19)  through  by  B  and  obtain  approximately 


vy 

^  •  <op  -r>  =  ^ 


a  (V  x  z) 
P  f 


which,  for  an  irrotational  neutral  wind,  becomes 
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and  where  G  is  the  combined  acceleration  due  to  gravity  and  the  neutral  wind,  as  defined 
after  Eq.  (8),  K  is  Boltzman's  constant,  and  Tg  and  are  the  electron  and  ion  temper¬ 

atures,  respectively. 

It  is  interesting  to  note  that  the  dynamo  current  appearing  in  Eq.  (40)  contains  a  sub¬ 
stantial  diffusion  contribution  (a  "ring  current"  in  the  case  of  axially  symmetric  ir¬ 
regularities)  but  that  its  divergence  is  zero.  Thus,  if  we  take  the  collision  frequency  and 

ion  mass  to  be  uniform  over  distances  of  interest,  so  that  Vo  <•>  Vn,  we  obtain  no  contri- 

P 

bution  to  the  final  cross  product  in  Eq.  (42)  from  the  diffusion  terra  in  V^.  We  are  left, 
then,  with 


(”p  “T")  •  ('A  -  f ) 


(43) 


Equation  (43)  is  close  to  the  final  form  of  one  of  the  two  main  MHD  equations  that  constitute 
the  Hele-Shaw  analogy,  as  we  shall  soon  see.  Before  turning  to  the  cell-fluid  equations, 
however,  we  must  complete  our  consideration  of  the  plasma  velocity.  We  return,  therefore, 
to  Eq.  (38),  which  may  be  rewritten  as 


Clearly,  for  conditions  of  interest,  the  cross-field  drift  produced  by  the  slow  gravi¬ 
tational  fall  and  even  by  the  neutral  wind  (in  this  low-collision-frequency  regime)  are 
•  •  •  * 

negligible  compared  with  the  E  x  B  drift.  More  importantly  for  our  present  purposes,  normal 
cross-field  diffusion  may  be  ignored  for  gradient  lengths  down  to  100  meters,  but  not  for 
very  sharp  gradients. 

It  is  interesting  and  pertinent  to  consider  in  detail  the  source  of  the  diffusion  term 

in  Eq.  (44).  Appearing  in  the  effective  electrostatic  potential  gradient  in  the  first  term 

of  Eq.  (38)  is  a  term  associated  with  ambipolar  diffusion  and  arising  from  the  pressure 

gradient  of  the  electron  gas.  This  term  also  is  included  in  the  left-hand  side  of  Eq.  (43). 

•  •  •  . 

The  relationship  between  this  term  and  the  imposed  electric  field,  E,  must  depend  on 
boundary  conditions  not  considered  here.  A  solution  for  E  in  terms  of  the  boundary  con¬ 
ditions  and  the  ambipolar  diffusion  force,  therefore,  is  relevant  to  the  assessment  of 
limitations  on  the  Hele-Shaw  analogue  as  well  as  for  a  full  understanding  of  striation 
dynamics. 

The  relevant  point  for  our  present  purpose  is  that  in  the  heavy-ion,  low-collision 
approximation  appropriate  for  the  simplest  Hele-Shaw  analogue,  the  electron-diffusion  con¬ 
tributions  from  the  two  terms  in  Eq.  (43)  cancel,  leaving  an  ion-diffusion  contribution  to 
the  cross-field  drift  described  by  Eq.  (44).  For  an  equilibrium  plasma  (T^  =  Tg),  the  upshot 
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is  diffusion  terms  of  equal  magnitude  but  opposite  sign  in  the  two  equations  we  need  for  the 
analogy.  The  diffusion  contribution  to  the  conduction  current  contained  in  parentheses  on 
the  left-hand  side  of  Eq.  (43),  then,  may  be  ignored  under  the  same  conditions  that  the 
diffusive  drift  in  Eq.  (44)  may  be. 

At  this  juncture,  then,  we  find  that  we  can  proceed  with  the  Hele-Shaw  analogue  without 
concern  about  diffusion  down  to  (night-time,  F-layer)  plasma  scales  on  the  order  of  100  m. 
We  note  that  this  restriction  is  quite  acceptable  for  pursuing  questions  of  outer-scale 
evolution,  but  the  general  question  of  diffusion  effects  must  be  regarded  as  open.  Mean¬ 
while,  we  turn  our  attention  to  dynamics  of  the  Hele-Shaw-cell  fluid. 

#  ,  ,  , 

We  note  that  a  plasma  undergoing  an  E  x  B  drift  in  an  approximately  uniform  magnetic 
field  moves  essentially  as  an  incompressible  fluid.  The  cell  fluid  is  incompressible,  and 
since  it  is  undergoing  two-dimensional  flow,  its  velocity  may  be  expressed  as  the  curl  of  a 
vector  potential  (Batchelor,  1970)  whose  only  non-zero  component  is  a  stream  function, 
^(x,y)  as  follows: 

u  =  V  x  (i|> z)  (45 ) 

or  u  =  V  ip  x  z.  (46) 

But  the  steady-flow  approximation  also  holds  in  the  cell,  so  we  can  solve  the  Navier-Stokes 
expression,  in  the  form  of  Eq.  (22b),  for  the  velocity  and  multiply  through  by  the  density  to 
obtain 


(47) 

(48) 

(49) 

(50) 


where  we  have  ignored  effects  on  the  observed  gap-center  velocity  due  to  the  cross-gap 
gradients  of  fluid  density  and  velocity.  That  the  velocity  is  not  uniform  across  the  gap  is 
stated  by  Eq.  (18),  and  it  follows  that  the  cross-gap  density  is  not  uniform  in  the  fluid- 


interchange  region.  The  effects  of  cross-gap  gradients  on  the  dynamics  of  u^  on  x-y  scales 
significantly  larger  than  d  are  thought  to  be  innocuous,  but  this  conjecture  has  not  been 
proven  theoretically  or  experimentally. 


C.  SUMMARY  OF  THEORETICAL  RESULTS 

The  first  point  of  analogy  between  magnetoplasmas  and  Hele-Shaw  fluids  is  the  similar¬ 
ity  between  the  plasma  momentum  equation  and  the  Navier-Stokes  equation  as  represented  by 
Eqs.  (22a  and  b).  The  main  caveat  on  the  latter  is  that  its  drag  term  depends  upon  an 
assumption  of  zero  density  gradient  across  the  gap.  Related  assumptions  arise  in  completing 
the  Hele-Shaw  analogue  of  the  magnetoplasma  electrodynamics.  These  and  other  pertinent 
assumptions  and  approximations  are  identified  in  Table  1. 


Table  1.  Pertinent  Assumptions  and  Approximations 


Hele-Shaw  Fluid 
Incompressible 
Two-dimensional 
No  surface  tension 
No  capillary  action 
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We  note  that  the  assumption  of  negligible  diffusion  in  the  magnetoplasma  restricts 
application  of  Hele-Shaw  analysis,  as  presently  understood,  to  plasma  scales  larger  than 
about  100  m.  This  is  quite  acceptable  for  outer-scale  considerations,  which  are  important 
for  assessment  and  prediction  of  radiowave  scintillation  effects. 

We  now  rewrite  the  central  results,  Eq.  (43),  (44),  (45),  and  (50),  in  light  of  the 
approximations  summarized  in  Table  1  and  invoke  the  continuity  equations  for  both  the 
magnetoplasma  and  the  cell  fluid  to  obtain  the  results  in  Table  2. 


Table  2.  Mathematical  Analogy 


He le- Shaw  Fluid 


Magnetoplasma 


u  =  V  ip  x  z 


*  (pVJ>)  =  ( 


* -ST/ *  2 


V  6  -+ 

M  -  B 

U  =  -  x  — 

B  B 


(-Y)= 


—  +  •  Vp  =  0 


Jl  +  •  Vn  =  0 


We  see  from  the  six  equations  in  Table  2  that  there  is  a  consistent  analogy  between  (1) 

the  local  electrostatic  potential,  scaled  by  the  magnetic  field  strength,  and  the  cell-fluid 

stream  function  (alternatively  between  the  electric  field  and  the  gradient  of  the  stream 

function);  (2)  the  magnetic  field  direction  and  the  cross-cell  gap  direction;  (3)  the  plasma 

drift  velocity  and  the  fluid  flow  velocity;  and  (A)  the  electron  number  density  and  the 

fluid  mass  density.  The  key  analogy  (A)  arises  from  one  between  the  Peterson  conductivity 

and  the  fluid  density,  given  uniform  collision  frequency  and  ion  mass  (as  well  as  magnetic 

G 

field  strength).  The  remaining  analogous  quantities  are  the  effective  dynamo  velocity,  — 

J2->  vn’ 

•  •  d  -* 

and  the  fluid  free-fall  velocity,  — ~  .  Gravity,  g,  in  the  laboratory  plays  the  role  of 

*  -* 

wind-augmented  gravity,  G,  in  the  field,  and  the  analogy  is  completed  by  the  drag  imposed 

8  r 

through  the  scaled  kinematic  viscosity,  — ,  in  the  cell  and  plasma-neutral  collisions,  v^, 
in  the  field. 

The  analogous  parameters  are  summarized  in  Table  3. 


Table  3.  Analagous  Parameters 


He  1  e-Shaw  Fluid 


Velocity,  uA 
Stream  function,  ip 


Magnetoplasma 

Velocity, 

Scaled  electrostatic 

d) 

potential,  -*■ 

D 

Field-aligned  unit  vector, 


Cross-gap  unit  vector,  z 


CD|0S+ 


III.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


A.  THE  CELL 

The  mathematical  analogy  between  fluid  in  a  Hele-Shaw  cell  and  a  dvnamo-convec ted 
magnetoplasma  is  enticing.  It  offers  an  opportunity  to  study,  by  means  of  a  conceptually 
simple  laboratory  device,  very  complicated  dynamic  processes  otherwise  accessible  onlv  by 
means  of  expensive  field  experiments  and  ambitious  computations.  In  practice,  however,  the 
experimental  procedures  are  not  without  their  own  problems.  Certain  practical  difficulties 
had  been  identified  in  early  experiments  by  Thomson  (1971)  and  by  Thomson  and  Hoc'nstiir. 
(1971).  Some  difficulties  were  alleviated  in  a  follow-up  set  of  experiments  by  Ozarski , 
Piccirelli,  and  Thomson  (1973). 

There  remained,  however,  several  needs  for  quantitatively  useful  experiments.  The 
early  experiments  had  suffered  from  an  inability  to  maintain  a  constant  gap  dimension,  d, 
across  the  face  of  the  cell,  which  resulted  in  nonuniforn  dynamics  via  the  free-fall  driving 

d  ■  g  i 

velocity,  -  nl  .  The  follow-up  experiments  had  attempted  to  solve  the  problem  bv  means  of  a 
on 

small  cell,  in  which  the  hydrostatic  head  did  not  cause  the  cell  to  bulge  appreciably  in  the 
middle.  The  smallness  of  the  cell,  however,  limited  the  horizontal  window  over  which  data 
could  be  analyzed  and  the  duration  of  experimental  runs  permitted  by  the  vertical  dimension. 
Moreover,  repeatability  of  fluid  behavior  in  the  early  cells  had  never  become  truly  satis- 
f ac  torv . 

In  an  attempt  to  realize  the  promise  of  the  Hele-Shaw  technique,  a  new  cell  and  more 
controlled  laboratory  procedures  were  developed  under  this  contract.  A  larger  cell  was 
designed  and  built,  employing  thicker  glass  for  stiffness  and  mechanically  supported  for 
convenient  and  repeatable  initiation  of  experimental  runs;  it  is  illustrated  in  Figure  1. 
The  overall  assembly  consists  of  the  cell  itself,  a  frame,  a  light  box,  and  a  floor  stand. 
The  cell,  proper,  consists  of  two  sheets  of  3/4-in  float  glass,  cut  and  ground  to  overall 
dimensions  of  20  by  50  in.  Surface  flatness  is  better  than  one  can  measure  with  a  24-in 
straightedge  and  a  .0015-in  feeler  gauge  (smoother  than  38.1  ,,m  over  61.0  cm). 

The  cell  gap  was  set  at  less  than  a  millimeter  by  means  of  0.32-in  shim  stock,  1  in 
wide,  placed  around  the  perimeter  of  the  glass.  The  dimensions  of  the  f 1 ui d-con t ai n i ng 
portion  of  the  cell,  therefore,  are  45.7  cm  by  1.22  m  by  813  Mm  (18  by  48  by  .032  in). 
Deformation  under  hydrostatic  load  was  calculated  as  91.4  urn  (11%  of  the  gap  dimension)  at 
the  worst  point,  and  this  estimate  was  confirmed  approximately  hv  measurement. 
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The  cell,  in  its  frame,  is  supported  by  a  stand  that  permits  rotation  from  one  vertical 
orientation  (the  stable  one  in  which  the  lighter  of  the  two  fluids  "floats"  on  the  heavier), 
used  for  filling,  to  the  opposite  vertical  orientation  (in  which  the  heavier  fluid  overlies 
the  lighter),  used  for  experimental  runs.  A  mercury  switch  mounted  on  the  cell  frame 
initiates  a  minute  and  second  (to  the  nearest  tenth)  counter  when  the  cell  is  rotated. 

Backlighting  of  the  cell  fluids  is  provided  by  ten  40-watt  fluorescent  discharge  tubes, 
48  in  long,  mounted  in  a  flat-white  box  that  provides  an  overall  illumination  field  of  48  by 
80  in.  The  lights  are  mounted  horizontally,  with  the  plane  of  the  light  source  16  in  from 
the  interior  glass.  Two  sheets  of  1/8-in  translucent  plexiglass  (Rohm  and  Nass  No.  2447W), 
placed  10  and  2h  in  from  the  interior  glass,  are  used  to  diffuse  the  light. 

b.  LABORATORY  PROCEDURES 

Early  experimental  runs  with  the  cell  shown  in  Figure  1  revealed  that  control  of  the 
cell  gap  to  winthin  11%  by  means  of  thick  glass  was  not  sufficient  to  provide  satisfactorily 
uniform  dynamics.  A  "droop"  in  the  fluid  interface,  caused  by  glass  bulge,  produced  two 
undesirable  effects:  (1)  large-scale  convection  downward  near  the  middle  and  upward  near 
the  edges  of  the  cell,  and  (2)  more  rapid  onset  of  the  Rayleigh-Tay lor  instability  near  the 
cell  center.  The  residual  glass  bulge  was  compensated  for  by  means  of  an  empirically  developed 
procedure. 

We  found  that  release  of  some  fluid  from  one  end  of  the  cell,  after  filling,  with  the 
other  end  closed  to  the  atmosphere,  provided  a  means  for  controlling  the  bulge.  The  released 
amount  required  to  provide  a  straight  interface  is  surprisingly  large  and  quite  sensitive  to 
details  of  the  apparatus,  but  it  is  also  repeatable.  Initially,  it  was  found  that  releasing 
90  ml  after  completely  filling  the  565-ml  cell  produced  a  satisfactory  interface.  After 
very  minor  modifications  were  made  to  the  filling  apparatus  (which  had  the  effect  of 
changing  the  overall  hydrostatic  head),  release  of  130  ml  was  found  necessary  to  produce  a 
staight  interface.  The  latter  amount  has  been  found  to  be  consistently  satisfactory  in  the 
absence  of  any  modifications  to  the  apparatus.  Figure  2  shows  examples  of  a  Rayleigh-Tay lor 
unstable  interface  with  (right)  and  without  (left)  control  of  the  hydrostatic  head. 

Several  fluids  have  been  used  for  studying  hydrodynamic  instability  in  Hele-Shaw  cells. 
The  theoretical  work  of  Lord  Rayleigh,  around  the  turn  of  the  century,  on  gravitationally 
driven  instability,  and  that  of  Taylor  (1950),  including  an  algebraically  additive  acceler¬ 
ation,  led  to  experiments  by  Lewis  (1950),  who  employed  air/water,  air/glycerine,  and 
benzene/water  combinations.  Saffman  and  Taylor  (1958)  addressed  the  question  of  instability 
at  a  viscosity  interface,  by  means  of  further  experiments  employing  glycerine  and  air. 
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Experiments  also  have  been  performed  on  a  water/oil  interface  for  application  to  petroleum 
from  porous  ground.  In  experiments  similar  to  ours,  Wooding  (1969)  worked  with  a  diffusing 
interface  between  the  water  and  a  potassium  permanganate  solution. 

Thomson  (1970)  recognized  the  mathematical  analogy  between  the  dynamics  of  the 
Rayleigh-Taylor  instability  as  manifested  in  a  Hele-Shaw  cell  and  the  gradient-drift  in¬ 
stability  in  high-altitude  plasmas.  He  noted  also  the  need  to  maintain  the  kinematic 
viscosity  constant  across  the  two-fluid  interface,  in  order  to  simulate  the  simple  iono¬ 
spheric  situation  of  uniform  ion-neutral  collision  frequency.  (It  is  this  same  condition  - 
uniform  collision  frequency  -  that  is  violated  by  a  nonuniform  interplate  gap  in  the  cell.) 

For  our  purposes,  we  sought  two  miscible  fluids  with  an  appreciable  difference  in  den¬ 
sity,  P  ,  but  with  essentially  identical  kinematic  viscosity,  n  ,  where  n  =  p/v  and  p  is 
viscosity.  These  requirements  were  approximated  by  Ozarski,  Piccirelli,  and  Thomson  (1973) 
by  using  distilled  water  for  the  lighter  fluid  and  salt  water  for  the  heavier  fluid,  with 
glycerine  added  to  each  to  achieve  the  required  viscosities.  We  first  attempted  to  repeat 
their  conditions  approximately.  In  order  to  obtain  the  maximum  density  difference,  we  added 
sufficient  NaCl  to  fresh  water  to  obtain  a  saturated  salt  solution.  We  then  experimented 
with  various  combinations  of  the  solution  and  glycerine  and  of  dionized  H2O  and  glycerine. 

After  considerable  systematic  trial  and  error,  we  established  a  pair  of  fluids  that 
produced  satisfactory  kinematic  parameters  with  reasonable  consistency.  For  the  light  fluid, 
we  employ  approximately  64%  deionized  water  and  36%  glycerine;  the  heavy  fluid  is  approxi¬ 
mately  70%  saturated  salt  water  and  30%  glycerine.  The  kinematic  properties  are 
temperature-dependent,  so  the  fluids  are  measured,  readjusted,  and  measured  again  for  each 
experimental  run.  During  this  iterative  process,  the  heavy  fluid  is  dyed  to  make  it 
visually  and  photographically  distinguishable  from  the  light  fluid.  Typically,  a  difference 
in  density  of  about  10%  is  achieved  between  the  dyed,  heavy  fluid  and  the  clear,  light  fluid; 
the  difference  between  their  kinematic  viscosities  is  held  to  about  1%. 

The  fluid-dying  process  consists  of  adding  a  small  amount  of  NaOH  to  make  a  basic 
solution  and  then  introducing  a  pH  indicator,  thymol  blue.  The  most  difficult  experimental 
parameter  to  control  consistently  is  the  opacity  of  the  dyed  fluid.  The  opacity  to  white 
light  increases  with  addition  of  thymol  blue,  until  a  contaminant  precipitates  out  of  the 
fluid.  We  have  not  been  able  to  identify  fully  the  conditions  that  produce  precipitation, 
and  it  does  not  appear  always  to  occur  at  the  same  optical  opacity. 

Our  procedure  is  to  add  thymol  blue  until  a  sufficient  opacity  (subjectively  judged  on 
the  basis  of  experience)  is  reached  and  to  filter  out  precipitate  if  necessary.  Useful 
opacity  is  reached  with  a  mix  of  about  34  g  of  thymol  blue  to  2  liters  of  heavy  fluid.  The 
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filtering  process  alters  the  kinematic  properties  somewhat,  so  dyeing  is  an  integral  part  of 
the  iterative  preparation  of  the  heavy  fluid. 

Fluid  dying  is  important  because  relative  opacity  directly  affects  the  signal-to-noise 
ratio  achieved  in  photographic  recording.  We  employed  a  red  photographic  filter  to  enhance 
the  contrast  between  the  blue  and  white-lighted  fluids.  A  series  of  exposure  tests 
(Fremouw,  and  Miller,  1980)  resulted  in  selection  of  1/125  s  at  f2.8  with  a  #25  filter  and 
Tri-X  Pan  film  (ASA  400)  as  our  standard.  Early  cell  runs  were  developed  in  Kodak  D-76  for  8 
minutes  at  68°F.  Later,  development  was  pushed  to  10  minutes,  which  became  our  standard. 

Some  of  the  exposure-test  frames  were  of  the  cell  in  its  stable  orientation.  Micro¬ 
densitometer  scans  vertically  across  the  fluid  interface  in  such  a  frame  were  used  to  define 
an  effective  signal-to-noise  ratio 


S 

N 


(51) 


where  £  is  the  mean  value  of  film  density  in  the  light-fluid  portion  of  the  frame,  h  is  the 
mean  value  in  the  heavy-fluid  portion,  and  o^  is  the  standard  deviation  of  the  former.  An 
example  of  such  a  scan  is  shown  in  Figure  3  (top).  The  tested  frame  yielded  S/N  *  33  dB . 
Figure  3  (bottom)  contains  a  horizontal  scan  across  the  structured  portion  of  the  right-hand 
data  frame  in  Figure  2.  It  shows  that  S/N,  as  defined  in  Eq.  (51),  is  not  quite  achieved  in 
practice,  probably  because  of  fluid  path-length  effects  to  be  discussed  in  Section  IV  A. 
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Figure  3.  Densitometer  scans  across  stable  fluid  interface 
(top)  and  structured  interchange  region  (bottom) 


IV.  ANALYSIS  PROCEDURES 


A.  DENSITOMETRY 

Figure  4  shows  nine  data  frames  exposed  during  a  cell  run  performed  under  the  standard 
fluid  and  photographic  conditions  described  in  Section  III  B.  The  run,  1(2 1 ,  was  performed 
on  10  October  1980  and  lasted  for  just  over  two  minutes.  The  densities  of  the  heavy  and 
light  fluids  were  1.207  and  1.103  g/ml,  respectively,  and  the  respective  kinematic  viscosi¬ 
ties  were  2.46  and  2.50  centistokes. 

Data  frames  such  as  those  shown  are  scanned  with  a  microdensitometer  (Photometric  Data 
Systems  Series  1000),  the  output  of  which  is  presented  as  film  density.  In  principle,  the 
record  one  wants  for  analysis  is  one  of  fluid  density  as  a  function  of  x-y  position  on  the 
cell  face.  If  all  shades  of  gray  on  the  film  were  uniquely  related  to  density  levels  in  the 
fluid,  then  the  specific  relationship  between  them  would  be  pertinent.  In  practice,  how¬ 
ever,  we  are  dealing  mainly  with  two  density  levels,  and  we  believe  that  only  negligible 
mixing,  due  to  molecular  diffusion,  takes  place  on  microscopic  scales  during  the  couple  of 
minutes  duration  of  an  experimental  run. 

Intermediate  shades  of  gray  do  develop  during  a  run,  but  we  believe  them  to  result  more 
from  fluid  thickness  effects  than  from  fluid  density.  That  is,  as  fingers  of  the  dark  and 
light  fluid  interpenetrate,  a  coating  of  the  other  fluid  originally  occupying  that  portion 
of  the  cell  is  left  on  the  glass  surface.  Indeed,  one  would  expect  such  an  effect  in  view  of 
the  cross-gap  velocity  profile  described  in  Eq.  (18).  Thus,  the  tips  and  edges  of  fluid 
structures  tend  to  have  less  contrast  than  interior  portions  of  the  structures  having  the 
same  centerplane  density. 

In  view  of  the  foregoing  complication,  there  would  be  little  value  in  trying  to 
identify  the  quantitative  relationship  between  dye  density  and  mass  density  in  the  heavy 
fluid;  for  practical  purposes,  we  take  them  to  be  proportional.  That  is,  if  is  the  mass 
density  of  the  heavy  fluid  and  that  of  the  light  fluid  and 


AP  =  Ph  ~  P,  (52) 

we  suppose  the  dye  density  and,  therefore,  its  optical  extinction  coefficient,  to  be  propor¬ 
tional  to  Ap. 
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Accordingly,  the  emittance  of  an  elemental  surface  area  in  the  heavy-fluid  portion  of 
the  cell  is 


-kApd  sec  0 

L  =  L  e 
o 


(53) 


and  the  corresponding  photographic  exposure  is 

„  -k-5pd  sec  0 

E  =  E  e 
o 


(54) 


where  Lq  and  Eq  are  the  emittance  and  exposure,  respectively,  of  an  elemental  area  on  the 
light-fluid  portion  of  the  cell,  k  is  a  proportionality  constant,  d  is  the  cell  gap  thick¬ 
ness,  and  0  is  the  angle  between  the  light  ray  from  the  elemental  surface  area  and  the 
central  axis  of  the  camera  lens.  The  film  density,  D,  is  related  to  the  exposure,  E,  by 
properties  of  the  film,  a  and  Y,  as  follows: 


D=  a  +  Y  log  E. 


(55) 


From  Eqs.  (54)  and  (55),  we  have 

(56a) 
(56b) 


D  =  a  +  i  (log  E  -  k  .'c  d  sec  log  e) 
=  b  -  c  ".f  sec 


where  b  and  c  are  constants  given  by 
b  =  a  +  Y  log  E 

o 

and  c  =  0.434k. 

Inverting  Eq.  (28b),  we  get 

Ap  =  cos  o  (57) 

c 

where  constants  b  and  c  are  established  by  the  photographic  setup  and  the  fluid  dye, 
respectively. 

In  principle,  one  could  employ  Eq.  (57)  to  calculate  An  to  within  the  constant  c.  That 
is,  0  could  be  established  from  the  experiment  geometry,  and  b  could  be  evaluated  from  the 
measured  film  density  in  the  light-fluid  portion  of  a  frame.  The  densitometer  record  of 
film  density  then  could  be  presented  as  a  record  of  film  density  in  units  of  (minus)  c. 
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Indeed,  one  could  even  adjust  the  scale  factor  for  different  runs  to  account  for  the 
different  k  (and,  therefore,  c)  values  arising  from  the  modest  inconsistency  in  opacity 
achieved  in  the  dying  process.  The  inconsistency  then  would  affect  only  S/N  and  not  the 
signal  levels  themselves. 

Our  interest,  however,  is  in  behavior  of  the  power  spectrum  of  the  density  variations, 

and  not  in  the  actual  values  of  fluid  density.  Our  immediate  objectives  can  be  met  by 

Fourier  transforming  the  records  of  film  density  themselves,  and  this  is  what  we  have  done. 

★ 

Figure  5  shows  two  densitometer  scans  from  run  #21.  The  upper  scan  was  taken  in  the  y 
direction  across  the  interface  prior  to  initiation  of  the  Ray leigh-Tay lor  instability. 
Referring  to  Eq.  (56),  it  shows  that  the  constant  b  was  approximately  1.74  film-density 
units.  The  gradual  drop  in  film  density  in  the  heavy-fluid  portion  of  the  scan  is  attributed 
to  the  increasing  path  length  through  the  translucent,  dyed  fluid  (the  sec  0  trend). 

The  lower  scan  in  Figure  5  was  taken  from  right  to  left  across  the  upper  portion  of  the 
fluid-interchange  region  (the  region  in  which  fingers  of  light  fluid  are  protruding  into 
heavy  fluid)  in  the  third  frame  (44.4  s)  of  Figure  4.  The  zero  level  at  the  extreme  left  and 
right  demarks  the  non-i lluminated  cell  frame,  and  the  curved  baseline  level  linking  them 
clearly  shows  the  sec  0  effect.  It  represents  a  weak  trend  on  scales  considerably  larger 
than  that  of  the  structures  of  interest,  and  we  have  not  corrected  for  it.  Indeed,  its 
clarity  in  Figure  5  attests  to  uniformity  of  illumination,  which  was  a  design  objective  in 
the  new  cell.  (The  trend  variations  evident  in  the  top  scan  of  Figure  3  resulted  from 
reflections,  which  were  eliminated  later.) 

Note,  in  addition,  the  difference  in  film  density  between  the  light-fluid  portion  of 
the  upper  scan  in  Figure  5  and  the  light-fluid  finger  tips  in  the  lower  scan.  As  alluded  to 
earlier,  this  difference  is  thought  to  arise  from  a  difference  in  light-fluid  thickness;  it 
decreases  the  S/N  ratio  in  the  interchange  region  by  six  dB  or  so  compared  with  that  computed 
from  a  scan  across  the  initial  interface. 

The  film  employed  was  in  6  by  cm  format,  with  the  result  that  1  mm  on  the  cell 
transforms  to  43  urn  on  the  film.  Thus,  the  smallest  relevant  size  on  the  film,  corresponding 
to  the  d  =  813  urn  gap  size,  is  35  pm.  The  circular  slit  used  on  the  densitometer  has  a 
diameter  of  10.6  pm,  which  was  stepped  in  overlapping  6-pm  steps  in  the  x  direction  (corre¬ 
sponding  to  horizontal  on  the  cell).  Thus,  the  densitometer  parameters  were  closely  matched 
to  the  smallest  relevant  size,  providing  about  three  independent  samples  in  that  distance 
(six  sample  points  on  the  abscissa  of  Figure  5). 


♦Keep  in  mind  that  the  scans  are  across  negatives  of  the  prints  shown  in  Figure  2. 
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SAMPLE  POINTS  -  horizontal 


Figure  5.  Calibrated  densitometer  scans  across  pre- instability 

interface  (top)  and  upper  portion  of  structured  inter¬ 
change  region  (bottom)  in  Run  #21. 


The  film  was  similarly  matched,  the  response  of  its  modulation  transfer  function  being 
rated  as  80%  at  30  cyles/mm  (approximately  one  cycle  per  imaged  gap  dimension).  Overall, 
the  expected  resolution  was  just  slightly  coarser  than  the  scale  on  which  third-dimension 
dynamics  might  be  expected  theoretically  to  produce  extraneous  effects. 

Typically,  50  x-direction  densitometer  scans  were  made  on  a  data  frame,  yielding  a  scan 
every  300  pm  in  the  y  dimension.  Ultimately,  two-dimensional  statistics  (in  both  configur- 
aton  and  Fourier  space)  are  of  interest.  Analysis  to  date,  however,  has  been  restricted  to 
the  time  evolution  of  one-dimensional  statistics  in  the  x  direction.  In  particular,  our 
attention  has  been  focused  on  development  of  a  one-dimensional  outer  scale  near  the  original 
interface  between  the  two  fluids. 

Densitometer  scans  from  near  the  middle  of  the  interchange  region  have  been  inspected 
visually  to  select  the  one  showing  most  nearly  a  balance  between  heavy  and  light  fluid, 
which  is  then  defined  as  the  center  scan  for  that  frame.  The  center  scans  for  six  of  the 
nine  frames  shown  in  Figure  2  are  presented  in  Figure  6.  The  main  evolutionary  behavior 
illustrated  is  an  increase  in  size  of  the  visually  dominant  structures,  to  which  we  shall 
return  in  Section  V, 

Figure  6  also  illustrates  a  problem  that  has  been  a  source  of  some  difficulty  in  the 
present  Hele-Shaw  experiment:  namely,  intermittently  unreliable  operation  of  the  densito¬ 
meter.  When  operating  properly,  the  Photometric  Systems  1000  has  been  very  suitable  for  our 
purpose  in  terms  of  resolution,  accuracy,  and  other  characteristics.  Unfortunately,  how¬ 
ever,  the  instrument  is  old  and  not  fully  reliable.  It  is  located  at  the  Boeing  Company, 
from  whom  we  purchased  operating  time,  and  Boeing  chose  not  to  perform  major  maintenance  on 
it  because  of  its  age  and  low  use  rate. 

There  are  two  aspects  of  the  densitometer  problem:  a  slowly  changing  output  bias,  and 
variable  noise  behavior.  Both  effects  are  illustrated  in  Figure  7,  which  also  shows  our 
basis  for  coping  with  the  former.  The  figure  shows  two  scans  taken  across  a  photographic 
step  tab  at  different  times  in  a  single  densitometry  session.  Comparison  clearly  shows  that 
the  instrument  bias  and  noise-fluctuation  levels  changed  between  the  scans  shown.  We  have 
found,  fortunately,  that  the  bias  drift  is  slow  compared  with  the  time  taken  to  perform 
densitometry  on  a  single  data  frame.  We  have,  therefore,  been  able  to  calibrate  the  set  of 
scans  from  each  frame  by  means  of  step-tab  scans  obtained  between  processing  of  each  pair  of 
data  frames.  The  film  density  used  for  display  and  analysis,  may  be  reliably  compared  from 
frame  to  frame  throughout  a  run  because  of  this  calibration  procedure.  Moreover,  response 
variations  from  scan  to  scan  are  unimportant  for  analysis,  since  absolute  density  is  not  of 
concern  to  us. 

More  difficult  to  deal  with,  however  has  been  the  intermittent  noise  contained  on  the 
densitometer  records.  The  noise  is  enigmatic.  Some  scans,  such  as  the  upper  scan  in  Figure 
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Figure  7.  Two  densitometer  scans  across  a  photographic  step 
tab,  showing  different  bias  and  noise-fluctuation 
levels . 


7,  the  second  (44.4  s)  and  last  (85.2  s)  in  Figure  6,  and  both  scans  in  Figure  5,  are  free  of 
the  added  noise.  Others,  such  as  the  lower  scan  in  Figure  7  and  the  third  (54.6  s)  scan  in 
Figure  6,  are  inundated  essentially  throughout.  Still  others,  such  as  the  first  (34.2  s), 
fourth  (64.8  s),  and  fifth  (75.0  s)  in  Figure  6,  are  partly  disturbed,  sometimes  in  the  same 
portion  of  the  scan  (leftmost  portion  of  the  first  and  fourth). 

The  scans  in  Figure  6  are  from  our  most  recent  densitometer  session,  and  its  shows  the 
most  seriously  degraded  data  in  our  collection;  the  problem  is  getting  worse.  Fortunately, 
the  spectral  characteristics  of  the  noise  experienced  have  kept  it  from  significantly 
affecting  our  outer-scale  investigation. 

B.  SPECTRAL  ANALYSIS 

Spectra  of  the  six  scans  appearing  in  Figure  6  are  shown  in  Figure  8.  Our  spectral 
analysis  is  performed  over  a  processing  window  of  2028  sample  points  approximately  centered 
on  the  cell  to  avoid  edge  effects.  The  spectra  in  Figure  8  were  derived  from  a  window 
spanning  sample  points  1000  through  3048  in  Figure  6.  Inspection  of  Figure  6  shows  that  the 
processing  window  was  essentially  free  of  the  spurious  noise  in  the  scans  from  44.4  s,  64.8 
s,  and  85.2  s,  while  the  54.6  s  scan  was  contaminated  throughout.  Using  the  44.4  s  scan  as  a 
noise-free  standard,  we  see  that  the  main  effect  of  the  contamination  on  the  partially 
affected  scans  (34.2  s  and  75.0  s)  is  to  raise  the  noise  floor  at  frequencies  above  about  120 
cycles/data  window.  The  effect  on  the  fully  contaminated  spectrum  (54.6  s)  is  similar  in 
character,  but  stronger,  extending  down  to  perhaps  85  cycles/window.  In  our  outer-scale 
analysis,  we  have  made  use  of  data  only  from  below  the  frequency  at  which  noise  appears  to 
dominate  each  spectrum. 

In  conjunction  with  our  analysis,  we  investigated  the  effect  of  various  processing 
procedures  on  data  obeying  a  cutoff  power-law  spectrum  (Fremouw,  Lansinger,  and  Miller, 
1979).  For  the  investigation,  we  generated  random  time  series  whose  power  spectral  density 
(PSD)  was  controlled  to  be  as  follows: 

P(K)  =  C  [l  +  (LK)2]  (58) 

in  which  the  low-frequency  power  spectral  density  (PSD),  C,  the  outer  scale,  L,  and  the 
high-frequency  spectral  index,  n,  could  be  controlled.  The  generation  routine  also  per¬ 
mitted  selection  of  Fourier  components  that  were  either  synchronous  or  asynchronous  with  the 
processing  window  and  whose  phases  were  either  ordered  or  random. 

Our  procedure  was  to  generate  a  time  series  with  selected  spectral  characteristics  and 
then  to  process  it  with  the  same  FFT  routine  used  for  analyzing  our  experimental  Hele-Shaw 
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data.  The  purpose  was  to  assess  the  relative  efficacy  of  specific  analysis  procedures,  such 
as  windowing  with  different  time-domain  tapers,  when  used  to  process  data  for  outer-scale 
studies. 

One  unexpected,  although  retrospectively  understandable,  behavior  was  noted  in  prepar¬ 
ation  for  our  tests.  We  found  that  spectral  leakage  depends  very  much  on  the  synchronous/ 
asynchronous  nature  of  spectra,  quite  independently  of  the  phase  relationship  between  the 
Fourier  components  contained.  Synchronous  data  whose  spectral  lines  lie  at  even  half-bins 
in  Fourier  space  preclude  spectral  leakage  even  for  a  rectangular  temporal  window  because 
the  corresponding  spectral  window  has  nulls  at  such  points.  Spectral  leakage  through  window 
sidelobes  is  maximized  by  placing  the  components  at  odd  half-bin  frequencies,  and  we  used 
such  tailored  synchronous  spectra  for  our  window  tests. 

For  our  tests,  we  spectrally  analyzed  a  data  set  with  the  foregoing  spectral  character¬ 
istics  eight  times  with  different  windows.  In  addition  to  rectangular,  the  windows  employed 
(Harris,  1978)  were  the  following:  B lackman-Harr i s  with  maximum  sidelobes  down  by  61,  67, 
74,  and  92  dB ;  and  cosine  with  10%,  20%,  and  30%  taper.  We  then  performed  three-parameter 
best  fits  between  the  resulting  numerical  spectra  and  the  analytical  spectrum  given  in  Eq. 
(58).  Our  concern  was  the  effect  of  data  processing  on  ascertaining  spectral  parameters  in 
real  data,  particularly  the  outer  scale.  To  test  for  biases  introduced  in  the  curve-fitting 
procedure,  we  processed  a  time  series  generated  from  synchronous  data  tailored  for  zero, 
rather  than  for  maximum  spectral  leakage. 

From  our  test  results,  we  were  able  to  draw  three  conclusions  about  processing  data 
that  obey  a  cutoff  power-law  spectrum.  First,  choice  of  an  FFT  window  and  of  whether  to 
perform  spectral  smoothing  is  not  of  great  consequence  in  identifying  relevant  spectral 
parameters,  but  there  are  some  differences.  The  best  efficacy  for  evaluation  of  the  outer 
scale  and  spectral  index  was  obtained  by  means  of  a  30%  cosine  taper,  and  we  subsequently 
used  such  a  window  in  our  processing.  Second,  mean  square  error  is  not  a  very  useful 
indicator  of  fit-parameter  efficacy. 

Third,  the  most  important  effects  revealed  in  our  tests  stemmed  from  the  curve-fit 
procedure,  and  they  can  be  rendered  innocuous  by  insuring  sufficient  iterations.  If  in¬ 
sufficient  iterations  were  employed  in  our  fit  procedure,  the  following  consistent  biases 
resulted:  spectral  strength,  C,  was  too  large;  spectral  index,  n,  was  too  small  (shallow); 
and  outer  scale,  L,  was  too  large.  The  foregoing  biases  are  all  understandable  for  three- 
parameter  fits  employing  data  from  the  entire  spectrum,  at  least  when  one  fits  in  loga¬ 
rithmic  coordinates.  A  three-parameter  fit  routine  suitable  for  use  with  Eq.  (58)  and 
providing  sufficient  iterations  in  reasonable  time  has  been  developed  by  McDonald  (private 
communication,  1979),  and  we  employed  it  for  our  analysis. 
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V.  RESULTS 


In  addition  to  showing  experimental  PSD's  obtained  from  Hele-Shaw  run  #21  by  means  of 
the  spectral  analysis  described  in  Section  IV  B,  Figure  8  shows  the  best  fits  thereto  of  Eq. 
(58).  As  attested  by  the  smooth  curves,  the  cutoff  power-law  form  provides  a  very  effica¬ 
cious  representation  of  the  experimental  spectra.  The  most  striking  feature  of  the  experi¬ 
mental  PSD’s  is  a  consistent  transition  from  a  low-frequency  white-noise  regime  to  a  high- 
frequency  power-law  tail.  We  identify  the  spectral  break  with  an  outer  scale  and  note  that 
it  invariably  occurs  in  Hele-Shaw  spectra. 

Recall  that  Run  #21  employed  laboratory  fluids  carefully  controlled  to  permit  relevance 
of  the  theoretical  development  summarized  in  Section  II.  Figure  9  shows  center-scan  spectra 
and  best  fits  for  a  run,  #7,  in  which  the  fluids  were  not  so  controlled.  The  densities  of 
the  heavy  and  light  fluids  were  1.20  and  1.09,  respectively,  so  the  difference  of  10%  was 
essentially  as  in  Run  #21.  While  the  kinematic  viscosities  in  Run  #21  were  very  nearly 
equal,  however,  those  in  Run  #7  were  very  different  (5.6  and  1.8  centistokes  for  the  heavy 
and  light  fluids,  respectively). 

Note  now  that  the  general  character  of  the  spectra,  including  the  existence  of  an 
outer-scale  break,  is  the  same  for  the  two  runs.  Essentially  the  same  instability  is 
responsible  for  structures  that  develop  at  the  two-fluid  interface  in  a  Hele-Shaw  cell 
whether  fluid  dissimilarity  is  dominated  by  density  (Wooding,  1969)  or  by  viscosity  (Saffman 
and  Taylor,  1958).  It  appears,  from  our  results,  that  the  general  character  of  the  struc¬ 
tural  spectrum  that  develops  is  quite  independent  of  the  specific  kinematic  properties  of 
the  fluids  employed. 

As  we  shall  soon  see,  the  dynamic  evolution  of  the  structure  and  its  spectrum  also  is 
qualitatively  independent  of  experimental  details,  but  not  quantitatively  so.  Before  pro¬ 
ceeding  to  a  discussion  of  spectral  evolution,  however,  we  present  one  more  representative 
set  of  center-line  scans  and  spectra. 

Figure  10  shows  densitometer  data  obtained  from  films  of  Run  #44  of  the  earlier  Hele- 
Shaw  experiments  performed  by  Thomson  (private  communication,  1979).  The  fluids  used  for 
this  run  were  similar  to  those  used  for  our  Run  #7,  in  that  the  kinematic  viscosities  of  the 
heavy  and  light  fluids  were  rather  dissimilar,  being  6.77  and  2.02  centistokes  respectively. 
The  density  difference  was  about  12%,  the  individual  densities  being  1.23  and  1.11  gm/1. 
The  cell  size  was  somewhat  different  from  that  described  in  Section  III  A.  The  fluid- 
containing  portion  measured  61  cm  by  46  cm  by  approximately  1.6  mm,  the  latter  (gap) 
dimension  varying  substantially  over  the  cell  face. 
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Thomson  employed  35-mm  film,  so  the  scaling  from  cell  to  film  is  different  from  that 
described  in  Section  IHb.  Moreover,  in  early  densitometry,  we  employed  a  46  inn  slit,  so  the 
resolution  is  coarser  than  in  the  runs  discussed  earlier.  (The  apparently  excellent  spec¬ 
tral  signal-to-noise  ratio  in  Figure  10  stems  from  a  fortuitous  matching  of  the  spectral 
rolloff  of  the  unresolved  noise  and  the  underlying  data;  the  spectra  are  not  meaningful 
beyond  about  200  cycles/data  window.) 

Another  difference  between  the  processing  of  Thomson's  Run  44  and  that  described  in 

Section  IV  is  that  detrending  had  to  be  included  in  the  former  to  remove  the  effects  of  non- 

* 

uniform  illumination  on  the  cell.  The  process  employed  was  developed  with  considerable  care 
(Fremouw  and  Miller,  1979),  but  any  detrending  raises  questions  regarding  spectral  behavior 
at  low  frequencies.  The  similarity  of  the  spectra  in  Figures  8  and  9,  for  which  detrending 
was  not  necessary,  to  those  in  Figure  10,  for  which  it  was,  lends  credence  to  the  latter.  In 
particular,  the  existance  of  a  spectrally  identifiable  outer  scale  is  very  clear. 

With  an  outer  scale  established,  we  can  move  to  investigation  of  its  evolution.  We 
have  employed  the  analytical  fits  shown  in  Figures  8,  9,  and  10  to  evaluate  the  center-line 
outer  scale  during  the  course  of  the  Hele-Shaw  runs  represented  in  them.  The  results  are 
presented  in  Figure  11  as  scatter  plots  of  L,  defined  in  Eq.  (58)  and  scaled  in  units  of 
distance  on  the  Hele-Shaw  cell  used  in  each  run,  as  a  function  of  time. 

For  the  fits,  we  have  employed  data  from  only  the  portion  of  the  experimental  spectrum 
judged  to  be  reliable  and  meaningful.  For  instance,  for  Run  #21,  we  have  used  only  the  bins 
between  1  cycle/data  window  and  the  frequency  at  which  spurious  noise  is  judged  to  start 
dominating  each  spectrum,  as  discussed  in  Section  IV  B.  For  the  noncontaminated  spectra 
from  Run  #21  and  for  the  other  data  sets,  frequencies  between  1  and  180  cycles/data  window 
were  used,  based  on  considerations  of  system  resolution  and  cell-gap  size. 

For  each  of  the  three  runs  represented  in  Figure  11,  we  have  superposed  a  linear  least- 
square  trend  line.  There  is  considerable  scatter  in  each  case,  but  the  sense  of  the  trend  is 
very  consistent;  outer  scale  increases  with  time  in  all  cases.  For  Run  #21,  some  of  the 
scatter  probably  results  from  a  residual  effect  of  the  spurious  densitometer  noise.  In 
particular,  the  third  data  point  (54.6  s)  may  be  erroneously  high  due  to  noise  effects.  Run 
#7  also  contains  an  outlying  point  (at  35  s).  In  this  case,  there  is  no  suspicion  of 
spurious  noise  contamination.  Rather,  the  departure  of  the  35-s  point  from  the  trend  stems 
from  statistical  inhomogeneity  in  a  finite  data  window.  Inspection  of  the  data  frame  and 

*It  is  detrending  that  produced  essentially  zero-centered  data  scans  in  Figure  10  as  con¬ 
trasted  with  Figure  6. 
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densitometer  scan  in  question  shows  that  the  spectrum  was  dominated  by  two  large  structures 
that  happened  to  fall  within  the  processing  window.  We  have  excluded  the  point  from  the 
trend  fit  for  Run  #7. 

The  key  point  to  note  in  Figure  11  is  that,  despite  the  point  scatter  for  both  experi¬ 
mental  and  statistical  reasons,  a  clear  and  consistent  trend  in  outer  scale  occurred.  In 
all  three  runs,  which  were  performed  under  rather  different  conditions,  an  outer  scale 
clearly  existed  and  it  evolved  from  smaller  to  larger  size  during  the  course  of  each 
experiment. 

On  the  basis  of  Hele-Shaw  data  analyzed  to  date,  one  could  only  speculate  about  the 
factors  controlling  the  growth  of  outer  scale.  We  note  that  for  the  two  runs  in  Figure  11 
for  which  kinematic  viscosity  was  not  held  constant  for  the  two  fluids,  #7  and  T44,  the 
growth  rate  obtained  from  the  best-fit  trends  in  outer  scale  are  rather  similar,  being  0.5 
and  0.4  mm/s  respectively.  A  considerably  larger  growth  rate,  1.4  mm/s,  was  obtained  when 
the  kinematic  viscosity  was  much  more  controlled  for  validity  of  the  theoretical  analogy,  in 
Run  #21. 
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VI.  CONCLUSION 


In  addition  to  providing  an  opportunity  for  relatively  inexpensive  quantitative  analy¬ 
sis,  the  Hele-Shaw  technique  permits  two-dimensional  visualization  of  structures  developed 
by  the  Ray leigh-Tay lor  instability.  One  may  watch  the  instability  develop  in  real  time,  in 
movie  format,  and  in  sequential  still  photographs. 

Except  for  a  few  early  runs  that  may  be  viewed  either  as  experimental  "disasters"  or  as 
"learning  experiences,"  we  have  found  that  the  qualitative  statistical  behavior  of  the  two- 
fluid  interchange  in  our  experiments  is  repeatable.  Moreover,  the  behavior  is  consistent 
with  that  experienced  in  the  more  successful  of  Thomson's  earlier  experiments.  The  quali¬ 
tative  development  may  be  described  in  terms  of  the  five  morphological  "stages"  identified 
in  Table  4. 


Table  4.  Observed  Morphology 

1.  First  Discernible  Structuring  —  earliest  notable  features,  usually  rising  voids; 

2.  Growth  of  Reciprocal  (Up  and  Down)  Fingers  With  Increasingly  Bulbous  Ends  and 

Filamentary  Stems  —  also,  development  of  lateral  warps; 

3.  Bifurcation  and/or  Trifurcation  of  Bulbous  Ends; 

4.  Complicated  Growth,  including: 

A.  Poly furcation  on  stem  sides, 

B.  Polyfurcation  on  bulbous  ends, 

C.  Cascading  (repeated  bifurcation  and/or  trifurcation  of  bulbous  ends); 

5.  Increasing  Occurrence  of  Quasi-Pinching  and  Merging. 

Table  4  is  ordered  chronologically,  but  the  stages  are  not  discrete.  Rather,  one  stage 
gives  way  to  the  next  with  no  clear  time  of  demarcation.  Stages  4  and  5  are  especially 
difficult  to  distinguish  from  one  another,  and  together  they  take  up  a  large  majority  of  the 
duration  of  a  run.  Stage  4-5  is  dominated  by  (1)  repeated  bifurcation  and  trifurcation  at 
the  tips  of  large-scale  fingers  of  penetrating  fluid  and  (2)  simultaneous  merging  of  old 
structures. 

The  aforementioned  processes  (1)  and  (2)  may  be  viewed  as  competing  for  control  of  the 
structural  spectrum.  Their  regions  of  dominance,  however,  are  somewhat  distinct.  Our 
quantitative  analysis  has  concentrated  on  the  region  near  the  original  fluid  interface. 
Merging  dominates  in  this  region,  as  documented  in  Section  IV.  The  bulbous  tips  of  the  large 
structures,  however,  are  continuously  spawning  new,  smaller  features.  In  these  regions,  one 
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would  probably  find  a  rather  different  one-dimensional  spectrum.  For  instance,  horizontal 
scans  near  the  bottom  of  the  heavy-fluid  penetration  region  reveal  clusters  of  small  finger- 
lets,  separated  by  relatively  large  voids  containing  only  low-density  fluid. 

If  the  foregoing  characterization  applies  to,  say,  the  equatorial  F  layer,  as  we 
believe,  then  it  is  clear  that  a  radio  wave  propagating  vertically  and  scanning  horizontally 
would  encounter  a  variety  of  local  regions  with  quantitatively  different  spectral  character¬ 
izations,  as  has  been  suggested  by  Wittwer  (private  communication,  1980).  It  may  be  that 
the  magnetoplasma  spectra  are  more  varied  than  the  cell-fluid  spectra.  Still,  the  simple 
characterization  in  terms  of  a  cutoff  power  law,  Eq.  (58),  that  we  have  found  so  consis¬ 
tently  efficacious  for  the  latter  can  provide  very  good  fits  to  ionospheric  spectra  obtained 
from  in-situ  probes,  as  illustrated  in  Figure  12. 

Figure  12  shows  an  in-situ  plasma  scan  and  the  corresponding  FFT  spectrum  obtained  from 
Atmospheric  Explorer  E  (Basu  and  Basu,  1978).  The  general  character  of  the  scan  is  very 
similar  to  that  of  Hele-Shaw  scans,  as  can  be  appreciated  by  comparison  with  the  most  noise- 
free  scan  (44.5  s)  of  Figure  6. 

In  the  upper  panel  of  Figure  12,  note  the  spectral  break  between  an  essentially  white- 

noise  regime  and  a  power-law  tail.  The  solid  line  superimposed  on  the  measured  spectrum 

represents  the  best  fit  by  Eq.  (58),  the  parameters  of  which  are  C  =  0.6  fractional  electron- 
.  .  —1  * 

density  units  per  km  ,  L  =  1.4  km,  and  n  =  2.9.  Other  scans,  such  as  that  shown  in  Figure  12, 
have  been  presented  by  Basu  et  al . ,  1979,  and  we  have  performed  fits  to  their  spectra.  In 
addition  to  the  spectrum  from  AE-E  orbit  7162,  shown  in  Figure  12,  we  performed  fits  to  two 
scans  from  orbit  7048,  two  from  orbit  7049,  and  one  from  orbit  7050.  The  best-fit  outer 
scales  obtained  were  1.3,  2.0,  1.6,  2.7,  and  1.4  km.  All  but  one  of  these  spectra  have  been 
published  recently  by  Basu  et  al .  (1980). 

Our  results  lead  us  to  conclude  that  (1)  the  cutoff  power  law  provides  a  very  suitable 
spectral  representation  of  both  He le-Shaw-f luid  and  magnetoplasma  data  and  (2)  an  effec¬ 
tively  finite  (smaller  than  experiment  data  windows)  outer  scale  can  be  defined  at  least 
locally  in  time  and  space.  Its  spatial  homogeneity  and  temporal  evolution  are  outstanding 
questions . 

We  submit  that  our  Hele-Shaw  experiments  have  shed  some  light  on  the  temporal  evolution 
of  the  outer  scale.  Our  results,  however,  are  limited  in  time  and  space  and  could  be 
extended  in  both  domains.  The  logical  first  step  in  that  extension  would  be  to  perform 

*We  have  not  assigned  any  significance  to  the  spectral  index,  n,  of  Hele-Shaw  data  because 
of  uncertainty  over  the  source  of  shades  of  gray  in  the  photographs.  See  Section  III  B. 
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experiments  in  the  existing  cell  with  a  different  (three-layer)  configuration  of  fluids. 

In  all  of  our  experiments,  the  outer  scale  has  continued  to  increase  until  the  runs 
have  been  terminated  by  fingers  reaching  the  top  and  bottom  of  the  cell.  Wittwer  (private 
communication,  1980)  has  suggested  that  the  continuous  increase  in  outer  scale  may  result 
from  the  effectively  infinite  reservoir  of  dissimilar  fluids  available  during  the  course  of 
a  run.  He  has  speculated  that  a  finite  reservoir  might  impose  a  limit  on  the  outer  scale.  A 
test  of  this  conjecture  should  be  possible  by  means  of  a  relatively  narrow  layer  of  heavy 
fluid  sandwiched  between  two  thicker  layers  of  light  fluid.  It  may  be  that  penetration  of 
the  heavy  layer  by  rising  fingers  of  light  fluid  would  signal  an  end  to  outer-scale  growth. 

In  the  closing  days  of  this  contract  effort,  we  attempted  to  produce  a  three-layer 
configuration  of  fluid  for  such  a  test.  We  encountered  experimental  difficulties  in  con¬ 
trolling  the  two  interfaces,  at  least  one  of  which  should  be  totally  nonstructured  at  the 
outset  of  the  run.  As  a  result,  the  runs  attempted  must  be  classified  as  "disasters" 
reminiscent  of  the  "learning  experiences"  encountered  at  the  beginning  of  our  two-layer 
experiments.  We  have  little  doubt  that  the  experimental  difficulties  could  be  overcome  by 
refinement  of  filling  procedures,  and  we  suggest  three-layer  experiments  as  a  potentially 
productive  endeavor. 

Whether  with  three-layer  or  resumed  two-layer  experiments,  a  second  endeavor  that  also 
could  prove  useful  would  be  to  explore  a  locally  homogeneous  characterization  of  the  spatial 
spectrum  in  regions  well  removed  from  the  initial  fluid  interface  (near  the  tips  of  heavy 
and  light  fingers).  Equally  productive,  we  suggest,  would  be  analytical  effort  directed  at 
quantitative  scaling  of  Hele-Shaw  results  to  magnetoplasma  dynamics,  together  with  experi¬ 
ments  designed  to  probe  the  control  of  outer-scale  growth  rate. 
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